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Abstract 

Background: Efficient utilization of both glucose and xylose is necessary for a competitive ethanol production from 
lignocellulosic materials. Although many advances have been made in the development of xylose-fermenting strains of 
Saccharomyces cerevisiae, the productivity remains much lower compared to glucose. Previous transcriptional analyses 
of recombinant xylose-fermenting strains have mainly focused on central carbon metabolism. Very little attention has 
been given to other fundamental cellular processes such as the folding of proteins. Analysis of previously measured 
transcript levels in a recombinant XR/XDH-strain showed a wide down-regulation of genes targeted by the unfolded 
protein response during xylose fermentation. Under anaerobic conditions the folding of proteins is directly connected 
with fumarate metabolism and requires two essential enzymes: FADH 2 -dependent fumarate reductase (FR) and Erolp. 
In this study we tested whether these enzymes impair the protein folding process causing the very slow growth of 
recombinant yeast strains on xylose under anaerobic conditions. 

Results: Four strains over-expressing the cytosolic (FRD1) or mitochondrial {OSM1) FR genes and EROl in different 
combinations were constructed. The growth and fermentation performance was evaluated in defined medium as well 
as in a complex medium containing glucose and xylose. Over-expression of FRD1, alone or in combination with EROl, 
did not have any significant effect on xylose fermentation in any medium used. Over-expression of OSM 1, on the other 
hand, led to a diversion of carbon from glycerol to acetate and a decrease in growth rate by 39% in defined medium 
and by 25% in complex medium. Combined over-expression of OSM 1 and EROl led to the same diversion of carbon 
from glycerol to acetate and had a stronger detrimental effect on the growth in complex medium. 

Conclusions: Increasing the activities of the FR enzymes and Erolp is not sufficient to increase the anaerobic growth 
on xylose. So additional components of the protein folding mechanism that were identified in transcription analysis of 
UPR related genes may also be limiting. This includes i) the transcription factor encoded by HAC1 ii) the activity of 
Pdilp and iii) the requirement of free FAD during anaerobic growth. 
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Background 

The yeast Saccharomyces cerevisiae is one of few yeast 
species capable of growing under strict anaerobic condi- 
tions [1]. This trait together with a high tolerance to- 
ward inhibitory compounds, e.g. ethanol, weak acids, 
phenolics and furaldehydes, has made S. cerevisiae the 
organism of choice for bioethanol production [2]. The 
development of sustainable processes for biofuel produc- 
tion is an important step in the efforts to reduce green- 
house gas emissions and becoming independent of fossil 
fuels [3]. The utilization of lignocellulosic raw material 
for the production of fuel-grade ethanol is one such 
process currently under development [4,5]. These raw 
materials are generated as waste in e.g. agricultural and 
forestry industries and contain a large fraction of fer- 
mentable sugars. The composition of the sugar fraction 
varies between different materials, but the largest part 
often consists of glucose and xylose [4]. Efficient utiliza- 
tion of both these sugars is necessary for the ethanol 
production process to be economically feasible [6,7]. S. 
cerevisiae is well-known for its capability to ferment 
hexoses, especially glucose. However, S. cerevisiae can- 
not naturally utilize xylose and thus has to be genetically 
modified. The simultaneous expression of the XYL1 and 
XYL2 genes from Scheffersomyces stipitis encoding xy- 
lose reductase (XR) and xylitol dehydrogenase (XDH), 
respectively, is one pathway that enables xylose utiliza- 
tion by S. cerevisiae [8]. Although many advances have 
been made in the development of xylose-fermenting S. 
cerevisiae strains, the efficiency is still far from matching 
that of glucose [9,10]. 

The ability of S. cerevisiae to grow under strict anaer- 
obic conditions is dependent on the activity of fumarate 
reductase (FR) enzymes [11,12]. S. cerevisiae has two 
genes, FRD1 and OSM1, encoding the cytosolic and 
mitochondrial FR enzymes, respectively [13,14]. These 
enzymes catalyse the reduction of fumarate to succinate 
while regenerating FAD from FADH 2 [15]. This is the 
last step of the reductive branch of the TCA-cycle where 
mitochondrial NAD is regenerated under anaerobic con- 
ditions (Figure 1) [16]. Single deletion of FRD1 or OSM1 
has no effect on the ability to grow under anaerobiosis, 
but a double deletion mutant cannot grow under such 
conditions unless an external electron acceptor is sup- 
plied (e.g. oxygen, menadione or phenazine methosul- 
fate) [11,12,17]. It has thus been proposed that FR 
enzymes provide the only way for yeast to regenerate the 
FAD/FMN prosthetic group of flavin enzymes that are 
required for growth under anoxia [11]. Erolp is one 
such essential flavin- containing oxidase, which normally 
uses oxygen as a final electron acceptor for FAD regen- 
eration. Erolp is a thiol oxidase [18,19] that operates to- 
gether with Pdilp [20] in the endoplasmic reticulum 



(ER) where they take part in the maturation of secretory 
proteins [21]. Both proteins are essential for the forma- 
tion of disulphide bonds and together they form a clas- 
sical proteinaceous electron relay system in which 
electrons are passed from the substrate polypeptide, via 
Pdilp and the membrane-bound Erolp, to eventually 
react with oxygen (Figure 1) [22-24]. 

Approximately one third of the yeast proteome is 
processed in the secretory pathway which makes this 
mechanism indispensable for cell viability, growth and 
function [21]. Many of the maturation processes taking 
place in the ER are regulated by a mechanism referred 
to as the unfolded protein response (UPR) under condi- 
tions of stress or increased folding requirements [25]. 
The ER contains a transmembrane protein encoded by 
IRE1 which has a cytosolic ribonuclease domain and a 
luminal sensor domain. The sensor detects rising levels 
of unfolded proteins inside the ER and activates the 
cytosolic domain. This activity removes an intron from 
HAC1 pre-mRNA and the fragments formed are ligated 
by Trllp ligase to form mature HAC1 mRNA. Haclp is 
a potent transcriptional activator of UPR target genes 
which include both PDI1 and EROl, and dozens of genes 
related to translocation, glycosylation, folding and deg- 
radation [26]. Thus, the UPR balances the biosynthetic 
rate of secretory proteins and fine tunes the ER capacity 
to meet cellular demands. 

Previous transcription analyses of recombinant xylose- 
fermenting strains have mainly focused on the changes 
that occur in genes encoding enzymes involved in cen- 
tral carbon metabolism [27-30]. Very little attention has 
been given to other fundamental cellular processes such 
as the folding of proteins. Recently we conducted a com- 
prehensive analysis of the dynamic changes in intracellu- 
lar metabolites that occur during batch fermentation of 
a glucose/xylose mixture using a recombinant XR/XDH- 
strain [31]. In that study we observed a marked accumu- 
lation of intracellular malate and fumarate during xylose 
fermentation. The direct connection between fumarate 
metabolism and protein folding under anaerobic condi- 
tions prompted us to further investigate this system dur- 
ing xylose fermentation. Due to the essential role of FR 
enzymes in anaerobic growth we hypothesized that the 
activity of FR is reduced on xylose which in turn would 
affect the capacity to regenerate the FAD required for 
proper folding of proteins in the ER. Such an impair- 
ment of the protein folding process could be the under- 
lying cause for the slow anaerobic growth on xylose. In 
this study we examined this hypothesis by over- 
expressing the two genes encoding FR, alone and in 
combination with EROl, and evaluated the effect on fer- 
mentation performance in two different media contain- 
ing glucose and xylose. 
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Figure 1 Metabolic pathways involving Frdlp, Osmlp, Pdilp and Erolp. Under anaerobic conditions theTCA cycle operates as two branches, 
a reductive path (red color) and an oxidative path (turquoise color). Osmlp catalyses the last step of the reductive branch of the TCA-cycle where 
fumarate is converted to succinate and FAD is regenerated from FADH 2 . Frdl p catalyses the same reaction in the cytosol. The regenerated FAD is used 
inside the ER during the maturation of secretory proteins. The Pdilp and Erolp enzymes mediate the oxidative folding of proteins in the ER. When a 
nascent polypeptide is translocated into the ER the reduced cysteines are oxidised by Pdi 1 p giving rise to a disulphide bridge in the protein. The 
reduced Pdilp is subsequently oxidised by the membrane bound Erolp which transfers the electrons to oxygen via covalently bound FAD. In vitro 
experiments have shown that free FAD can act as the electron acceptor in the absence of oxygen. 



Methods 

Strains and culture conditions 

Plasmids and yeast strains used in the study are listed in 
Table 1. Yeast strains were recovered from 20% glycerol 



stocks stored at -80°C on solid YPD plates (10 g/L yeast 
extract, 20 g/L peptone, 20 g/L glucose, 15 g/L agar) for 
two days at 30°C Yeast cultures were grown in liquid 
YPD medium for 14-16 h, or less when required, at 30°C 
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Table 1 Strains and plasmids used in this study 


Name 


Relevant genotype 


Reference 


Plasmids 






P^-ZD 1 tr 


/ IP A 3 


roil 
pzj 


ripiac 1 zo 


1 PI 11 




YipidCz I I 


1 IDA 3 


pij 


YIpDRI 


Ylplacl 28:TDH3p-GXF1-CYCh 


[34] 


YlpDR7 


Y[plac211:7DH3p-XV/.7(N272D)-/\DH7t; 
PGK1p-XYL2-PGK1t 


[35] 


YIpBBI 


Y\p\ac128:TDH3p-FRD1-CYCh 


This study 


YlpBB2 


Y\p\acU8:TDH3p-OSMl-CYClt 


This study 


YlpBB3 


Y\p\ac]28:TDH3p-FRD1-CYCh; TEFlp-EROl-EROlX 


This study 


YlpBB4 


Y\p\acl28:TDH3p-OSM1-CYCh; TEFlp-EROl-EROlX 


This study 


Yeast strains 




TMB3043 


CEN.PK2-1C; gre3-A; his3::PGK1p-XKS1-PGKh, HIS3; 
tal 1 ::PGK1 p-TAL 1 -PGK1X; tkl 1 ::PGK1 p-TKL 1-PGK1X; 
rki 1 ::PGK1 p-RKI 1 -PGKlt; rpe 1 ::PGK1p-RPE] -PGK1X; 
ura3, Ieu2 


[36] 


TMB3452 


TMB3043; ura3::YlpDR7; Ieu2 


This study 


TMB3455 


TMB3452;/eu2::Ylplac128 


This study 


TMB3456 


TMB3452; /eu2::YlpBB1 


This study 


TMB3457 


TMB3452; /eu2::YlpBB2 


This study 


TMB3458 


TMB3452; /eu2::YlpBB3 


This study 


TMB3459 


TMB3452; /eL/2::YlpBB4 


This study 



and 180 rpm in an orbital shaker. Transformants were 
selected on solid Yeast Nitrogen Base (YNB) medium 
(6.7 g/L YNB with ammonium sulphate and without 
amino acids) with 20 g/L glucose as carbon source. Leu- 
cine was supplemented when required at 220 mg/L. Cul- 
tivations in shake flasks had a culture volume of 10% or 
less of the flask volume. The defined medium used in 
fermentation experiments contained 2X YNB, 20 g/L 
glucose, 50 g/L xylose and 50 mM potassium hydrogen 
phthalate buffer at pH 5.5. The complex medium used 
in fermentation experiments contained 10 g/L yeast ex- 
tract, 20 g/L glucose, 50 g/L xylose, 50 mM potassium 
hydrogen phthalate buffer at pH 5.5, 10 g/L (NH 4 ) 2 S0 4 , 
1 g/L MgS0 4 and 6 g/L KH 2 P0 4 . All media used in an- 
aerobic cultivations were supplemented with Tween 80 
(400 mg/L) and ergosterol (10 mg/L). 

E. coli strain NEB5a (New England Biolabs, USA) 
was used for sub-cloning of plasmid DNA. Transfor- 
mants were selected on solid LB plates (5 g/L yeast ex- 
tract, 10 g/L tryptone, 10 g/L NaCl, 15 g/L agar, 
pH 7.0), supplemented with 100 mg/L of ampicillin, 
for 16 h at 37°C. Cultures of transformed E. coli were 
recovered from 25% glycerol stocks stored at -80°C 
and grown in liquid LB medium, supplemented with 
ampicillin, for 14-16 h at 37°C and 180 rpm in an or- 
bital shaker. 



Cloning procedures 

Standard molecular biology techniques were used for all 
cloning procedures [37]. Plasmid purification was per- 
formed using Genejet Plasmid Mini-prep Kit (Thermo 
Scientific, USA). DNA fragments digested with restric- 
tion enzymes were excised from agarose gel and purified 
using QIAquick Gel Extraction Kit (Qiagen, Germany) 
when needed. DNA fragments amplified by PCR were 
purified using Genejet PCR Purification Kit (Thermo 
Scientific, USA). All enzymes for PCR amplifications, li- 
gations and restriction digestions were obtained from 
Thermo Scientific unless stated otherwise. PCR products 
were amplified using the High Fidelity Phusion Hotstart 
II Polymerase or High Fidelity Taq Polymerase Mix 
(Thermo Scientific, USA). Competent yeast cells were 
prepared and transformed according to the PEG/LiAc 
method [38] and heat shock competent E. coli cells were 
prepared according to the Inoue method [37] and trans- 
formed according to the supplier s instructions. 

The gene FRD1 was amplified from yeast genomic 
DNA by PCR using primers FRDl-Fl-Spel (5-GAC 
TAG TAA ATG TCT CTC TCT CCC GTT G-3') and 
FRDl-Rl-Sall (5 -GTT GTC GAC TTA CTT GCG 
GTC ATT GGC-3). The gene OSM1 was amplified 
from yeast genomic DNA by PCR using primers OSM1- 
Fl-Spel (5-GAC TAG TAA AAT GAT TAG ATC TGT 
GAG-3') and OSMl-Rl-Sall (5 -GTT GTC GAC TCA 
GTA CAA TTT TGC TAT G-3 ). The amplified FRD1 
and OSM1 genes were digested with Spel and Sail and 
inserted between the TDH3 promoter and CYC1 termin- 
ator of the previously constructed vector YIpDRI, yield- 
ing the two plasmids YIpBBI and YIpBB2 (Table 1). The 
cloned genes were verified by DNA sequencing. 

The EROl gene, including a 300 bp downstream se- 
quence, was amplified from yeast genomic DNA by PCR 
using primers EROl-Fl-BamHI (5 '-GAG GAT CCA 
AAA TGA GAT TAA GAA CCG CCA TTG-3) and 
EROl-Rl-XhoI (5 -GTC TCG AGC CCT TGA AGA 
TGG TAC C-3'). The amplified EROl ORF and termin- 
ator was digested with BamHl and Xhol and inserted be- 
tween the TEF1 promoter and CYC1 terminator of the 
vector p426TEF (Table 1). The TEFlp-EROl-EROlt cas- 
sette was amplified from plasmid p426TEF-ER01 using 
primers p426TEF-Fl-SacI (5'-GGA GCT CAT AGC TTC 
AAA ATG-3') and p426ER01t-Rl-NdeI (5 '-CAT ATG 
CCC TTG AAG ATG GTA C-3'). The amplified cassette 
was digested with BamHl and Xhol and ligated into plas- 
mids YIpBBI and YIpBB2, digested with the same enzymes, 
yielding vectors YIpBB3 and YIpBB4, respectively (Table 1). 
The cloned cassette was verified by DNA sequencing. 

Strain construction 

S. cerevisiae strain TMB3452 was constructed by trans- 
forming strain TMB3043 with YIpDR7, linearized with 
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EcoRN. The integrative plasmids YIpBBl-4 were linear- 
ized with EcoRV, Aflll, Eco91I and AflLL, respectively, 
and used to transform strain TMB 3452, yielding strains 
TMB3456 (ERD1), TMB3457 (OSM1), TMB3458 (ERD1 
EROl) and TMB3459 (OSM1 EROl) (Table 1). The con- 
trol strain TMB3455 was constructed by transforming 
TMB3452 with YIplacl28, linearized with EcoRV. Cor- 
rect integration of the vectors was verified by PCR using 
primer pairs that annealed on the plasmid and on the 
chromosome outside the homologous regions of URA3 
or LEU2 loci. The functionality of the over-expressed 
FRD1 and OSM1 genes was verified by measuring the 
specific fumarate reductase activity in cell-free protein 
extracts from strains TMB3455-59 grown anaerobically 
on xylose. Strains TMB3456-59 had between 9- and 30- 
fold higher activities than the control showing that the 
over-expressed genes generated increased levels of func- 
tional enzymes (Additional file 1: Figure SI). The func- 
tionality of the over-expressed EROl gene was verified 
by determining the sensitivity toward DTT. Single col- 
onies of strains TMB3455-59 were inoculated in 5 mL 
YNB medium with 2% glucose and grown over-night. 
200 \iL of this culture was used to inoculate fresh 
medium and cultivated until the cells reached OD 620 nm ~ 1. 
A volume containing 2 x 10 6 cells (150-200 \iL) was 
spread on solid YPD medium after which a filter disc was 
placed in the centre of the plate and soaked with 20 \iL or 
30 (iL of 1 M DTT solution. The plates were incubated 
48 h at 30°C before measuring the diameter of the halo 
formed around the disc. The halos formed with strains 
TMB3458 and TMB3459 were 20-30% smaller compared 
with the other strains showing that the over-expressed 
EROl gene generated increased level of a functional en- 
zyme (Additional file 1: Figure S2). 

Anaerobic batch fermentations 

Yeast strains were pre-grown in YNB medium with 2% 
glucose until mid-exponential phase and inoculated into 
50 mL medium with 2% glucose and 5% xylose in 
250 mL crimp-sealed serum flasks at a concentration of 
0.04 g CDW/L. Prior to inoculation the headspace was 
flushed with N 2 -gas containing less than 5 ppm oxygen 
(AGA Gas AB, Sweden) for at least 30 min to ensure 
anaerobic conditions. Cultivations were carried out in 
a 30°C water bath with continuous stirring at 180 rpm 
using magnetic stirrers. To prevent diffusion of oxy- 
gen into the culture during sampling, the excess gas 
was released through a 0.2 \im filter connected to 
Norprene tubing submerged under water. Cell dry 
weight was calculated from OD measurements using a 
pre-determined relationship between OD and cell dry 
weight. Fermentation experiments were performed in 
biological duplicates. 



Preparation of cell-free protein extracts 

Yeast strains were cultivated anaerobically in crimp- 
sealed serum flasks as described above using 2X YNB 
medium with 2% glucose and 5% xylose. After 36 h, 
when xylose was the only carbon source, the cells were 
harvested by centrifugation at 4000 rpm and 4°C for 
10 min. The resulting cell pellets were resuspended in 
IX Lysis buffer (20 mM K,Na-P0 4 , 20 mM EDTA, 
1 mM DTT, pH 7.0) to obtain a final suspension of 1 g 
wet cells/mL. Protein extracts were prepared using a 
Precellys @ 24 homogenizer with cooling provided by the 
Cryolys® unit (Bertin Technologies, France) and Precellys 
Lysing Kit VK05 containing 0.5 mm glass beads. The 
unit was operated for 3 x 45 s at 5000 rpm with 30 s 
pauses between each cycle and at max 8°C. Cell debris 
was pelleted by centrifugation at 14,000 x g and 4°C for 
15 min and the supernatant was transferred to a small 
glass vial. Before sealing the vial the headspace was re- 
placed with N 2 -gas in an anaerobic chamber. The pro- 
tein extracts were stored at 4°C until the following day. 
Protein concentrations were determined with the Coo- 
massie Protein Assay Reagent (Thermo Scientific, USA) 
using BSA as standard (Thermo Scientific, USA). 

Fumarate reductase activity assay 

Fumarate reductase activity was determined in a 1.2 mL 
reaction mixture containing: 50 mM K,Na-P0 4 buffer 
(pH 7.3), 10 mM fumaric acid (pH 7.3), 0.4 mM FMN- 
Na, 2% (w/v) sodium hydrosulfite. All reagents, except 
protein extract, were added to a rubber-sealed cuvette 
inside an anaerobic chamber. The assay buffer was pre- 
pared by dissolving 2.762 g of KH 2 P0 4 and 2.837 g of 
Na 2 HP0 4 -2H 2 0 in 50 mL anaerobic water (separately) 
inside the anaerobic chamber, creating one 10X solution 
and one 2X solution, respectively. A 200 mM stock solu- 
tion at the correct pH was made by mixing 5 mL 10X 
with 25 mL 2X and completing the volume to 50 mL 
with anaerobic water. A 40 mM stock solution of fu- 
maric acid was prepared by dissolving 0.464 g fumaric 
acid in 50 mL of 5 mM K,Na-P0 4 buffer at pH 7.3. The 
pH was adjusted to 7.3 using 5 M NaOH and the vol- 
ume was completed to 100 mL with water. Solutions of 
the cofactor and reducing agent were prepared fresh be- 
fore each assay. Both compounds were weighed and 
brought into the anaerobic chamber. The FMN cofactor 
was dissolved in anaerobic water and sodium hydrosulfite 
was dissolved in 100 mM K,Na-P0 4 buffer, but only after 
all other reagents had been added to the cuvette to ensure 
that FMN becomes completely reduced to FMNH 2 . The 
reaction was started by addition of protein extract using a 
gas-tight syringe. The formation of FMN was measured at 
445 nm in cuvettes with light path length of 1 cm using 
an Ultrospec 2100 Pro spectrophotometer (Amersham 
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Biosciences Corp., USA) controlled by computer soft- 
ware SWIFT Reaction Kinetics v. 2.05 (Biochrom Ltd., 
Cambridge, UK). The absorbance was converted to con- 
centration units using the molar extinction coefficient 
of FMN (8 445 nm = 12,500 L/mol/cm). The slope was de- 
termined within an interval of at least 1 min, using at 
least three different dilutions of the protein extract. 
Only the slopes within the range 0.045-0.40 A/min were 
considered acceptable. One unit is defined as the con- 
version of 1 umol of fumarate per minute. 

Biomass determination and analysis of metabolites 

Optical density was measured at 620 nm using an 
Ultrospec 2100 Pro spectrophotometer (Amersham Bio- 
sciences Corp., USA). Cell dry weight was measured in 
triplicate by filtering a known volume of the culture 
through a pre-weighed nitrocellulose filter with 0.45 um 
pore size. The filters were washed with three volumes of 
water, dried in a microwave oven and weighed after 
equilibrating to room temperature in a desiccator. 

Concentrations of glucose, xylose, xylitol, glycerol, suc- 
cinate, acetate and ethanol were analysed by HPLC 
(Waters, USA). Aminex HPX-87H ion exchange column 
(Bio-Rad, USA) was used at 45°C with a mobile phase of 
5 mM H 2 S0 4 at a flow rate of 0.6 mL/min. All com- 
pounds were detected with a RID- 10 A refractive index 
detector (Shimadzu, Japan). 

Calculation of metabolic rates 

In defined medium (YNB) the volumetric production 
and consumption rates (q t and q s ), the maximum specific 
consumption rates of xylose (r xy j) and glucose (r g i c ), as 
well as the biomass yield per consumed substrate (Y x/s ) 
were calculated through linear regression according to 
Equations 1-4. The product yields per consumed sub- 
strate (Y i/s ) were calculated using the determined volu- 
metric rates according to Equation 5. The product yields 
per biomass (Y i/x ) were calculated as the ratio between 
Y i/s and Y x/s when xylose was the carbon source and ac- 
cording to Equation 6 when glucose was the carbon 
source. Due to the non-exponential growth on xylose, the 
maximum specific growth rate (^ max ) was calculated as 
the product between r xy i and Y x/s according to Equation 7. 
For comparative reasons // max on glucose was calculated 
in the same way. These values corresponded well with 
those calculated through regression of ln(OD) vs. time. 
The specific production rates were calculated as the prod- 
uct between // max and Y i/x according to Equation 8. 



r g ic 



dCi 
dt 



dC s 
dt 



V xyl '■ 



d(C X yl/ Cycj 

dt 



(1) 



(2) 



d(dCgi c /dt) 
dC x 

dC x 



Yi/s 
Yi/x 



dC S: 

Is 

dCi 
dC x 

= r.-Y. 



x/s 

^max ' Yi/x 



(3) 
(4) 
(5) 

(6) 

(7) 
(8) 



In complex medium all calculations were made as 
described above except for the maximum specific con- 
sumption rate of glucose which could not be determined 
accurately enough. Hence, the ^ max under these condi- 
tions was determined as the slope in a linear regression 
of ln(OD) vs. time and r g i c was calculated using Equation 7. 
The significance of a difference between two mean values 
was calculated using a two-tailed t-test assuming equal 
variance of the two samples with two degrees of freedom 
(m = n 2 = 2) unless stated otherwise. Balances of carbon 
and degree of reduction, determined for the glucose 
and xylose phases using calculated yield coefficients, are 
given as supplementary information in Additional file 2: 
Table SI. 

Transcriptional analysis of UPR related genes 

The list of genes reported as targets of the UPR were 
taken from [26]. Transcript data from cultivations of a 
XR/XDH-strain using glucose or xylose in aerobic and 
anaerobic conditions was taken from [27]. UPR target 
genes that had a significant change in expression (P< 
0.1) during the transition from aerobic to anaerobic con- 
ditions with xylose as carbon source were included in 
the analysis. Two genes commonly used as targets of the 
UPR (KAR2 and INOl) [27] did not pass the statistical 
requirement and were thus not included in the analysis 
(Additional file 2: Table S5). Hierarchical clustering was 
performed using the clustergram function in MATLAB 
R2010b with log 2 values of transcriptional changes. 
Euclidean distance was used to calculate the distance be- 
tween the data points and average linkage was used to 
generate the dendrogram. Data used for the hierarchical 
clustering is given as supplementary information in 
Additional file 2: Table S2. 

Results and discussion 

Genes targeted by the unfolded protein response are 
differentially expressed on xylose compared to glucose 

A previously performed transcripome study in a recom- 
binant XR/XDH-strain [27] was used to investigate the 
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expression level of genes associated with the UPR. It re- 
vealed that IRE1 and HAC1 levels were -65% and -50% 
higher (P = 0.02 and P = 0.05, respectively) in cells grow- 
ing anaerobically on xylose compared to cells growing 
aerobically (Figure 2). However, the transcript level of 
TRL1, encoding the ligase, was -50% lower (P = 0.02) 
(Figure 2), suggesting that HAC1 is mainly present as 
fragments rather than mature mRNA molecules. The 
transcript levels of these genes did not change signifi- 
cantly during the same transition (i.e. from aerobic to 
anaerobic conditions) on glucose (Figure 2) indicating 
that the response was associated with the type of carbon 
source available. The cluster analysis illustrated in 
Figure 3 indeed shows that a majority (61%) of the in- 
cluded UPR target genes were down- regulated on xylose 
when the conditions changed from aerobic to anaerobic 
(Group A). The genes belonging to this group represent 
functions taking place in all steps of the maturation 
process (translocation, glycosylation, disulphide bond 
formation, protein folding, vesicle budding and vesicle 
transport to the cell membrane) (Additional file 2: Table 
S2) and do not change during the same transition on 
glucose. The smallest group (Group B) contained those 
genes which are regulated equally on glucose and xylose, 
with only small differences in abundance between the 
two carbon sources. A small number of genes were up- 
regulated on xylose while remaining largely unchanged 
on glucose (Group C). These genes are related to peroxi- 
some biogenesis (PEX4), ubiquitin recycling (DOA4) and 
the cell integrity signalling pathway (TUS1). The three 




-2 J 1 

IRE1 HAC1 TRL1 

Figure 2 Transcriptional response of genes involved in the 

regulation of the UPR. Bars show the log 2 fold change of the 

indicated genes during a transition from aerobic to anaerobic 

conditions in a XR/XDH-strain utilizing xylose (red bars) or glucose 

(turquoise bars) as carbon source [27]. Values are given as mean fold 

change of two independent measurements in each conditions and 

errors represent the 95% confidence interval of the mean {df=2). 
\ J 



genes ECM3, ECM8 and YET2 encode proteins with un- 
known functions. RIB1 encodes GTP cyclohydrolase II 
which catalyses the first step of the riboflavin (FAD and 
FMN) biosynthesis pathway. The fourth group contains 
genes which are strongly up-regulated in response to an- 
aerobiosis in both conditions (Group D). However, in 
contrast to the genes in Group B, the transcript levels of 
these genes were in general lower on xylose compared 
to glucose. This analysis shows that the UPR responded 
differently to xylose compared to glucose. Furthermore, 
the repression of a majority of the target genes suggests 
that the capacity of the secretory pathway could be lim- 
ited during xylose fermentation. 

Consistent with previously reported expression pat- 
terns on non-fermentable carbon sources and entry into 
stationary phase [39,40], the expression of EROl was 
reduced by -30% (P < 0.02) on xylose compared to glu- 
cose, both under anaerobic and aerobic conditions 
(Additional file 2: Table S5) [27]. Hence, a reduced cap- 
acity to form disulphide bonds in nascent polypeptides 
translocated to the ER could be one underlying cause for 
the low anaerobic growth rate on xylose. The previously 
reported accumulation of intracellular fumarate during 
xylose fermentation [31] suggested that recombinant 
XR/XDH strains have a limited FR activity under such 
conditions, which would lead to insufficient recycling of 
the FAD co-factor required by Erolp. This hypothesis 
was investigated by evaluating the physiological effect of 
over-expressing the cytosolic and mitochondrial FR en- 
zymes, alone and in combination with EROl, during fer- 
mentation of a mixture of glucose and xylose. 

Over-expression of cytosolic FR has limited impact on 
xylose fermentation 

Due to the higher contribution of the cytosolic FR to the 
total cellular reductase activity the over-expression of 
FRD1, alone and in combination with EROl, was evalu- 
ated first. When cultivated in defined medium (2X YNB) 
with 20 g/L glucose and 50 g/L xylose, TMB3456 (FRD1) 
and TMB3458 (FRD1 EROl) consumed the glucose within 
20 h and ca. 70% of the xylose within 90 h (Figure 4B 
and D). The two strains behaved very similar to the con- 
trol strain (Figure 4 A) and no significant physiological 
differences from this strain were observed (Additional 
file 2: Table S3). These results suggest that the imple- 
mented genetic modifications had no positive effect on 
the anaerobic growth on xylose in defined medium. 
However, the modifications only ensured that the genes 
required for regeneration of FAD co-factors were present 
(Figure 1). If substrates are not available in sufficient 
amounts (in this case unfolded polypeptides) these modi- 
fications are not effective in regenerating the co-factor. 
The poor repressive capability of xylose leads to an up- 
regulation of MDH2 [27] (Additional file 2: Table S5), 
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Figure 3 Hierarchical clustering of genes induced by the UPR. Genes that are induced by the UPR [26] and changed significantly (P< 0.1) 
during a transition from aerobic to anaerobic conditions in a XR/XDH-strain utilizing xylose [27] were analysed by hierarchical clustering. The heat 
map was generated using the log 2 fold change between the different conditions. The analysis shows that a majority (61%) of the included UPR 
target genes are down-regulated on xylose when the conditions change from aerobic to anaerobic (Group A). These genes do not change during 
the same transition on glucose indicating that the UPR responds differently to xylose compared to glucose. 



encoding cytosolic malate dehydrogenase, which has pre- 
viously been suggested as an unbeneficial reaction that 
directs carbon away from amino acid synthesis by con- 
verting oxaloacetate into malate [31]. A reduced synthe- 
sis of amino acids could thus hamper the formation of 
polypeptides and mask the effect of the genetic modifica- 
tions. To investigate whether the strains were limited by 
amino acid synthesis they were cultivated in a complex 
medium containing yeast extract. Changing from a de- 
fined medium to a complex medium led to several sig- 
nificant changes in the physiology of all strains, e.g. 
increased maximum growth rate, increased biomass yield 
and reduced by-product yields (Additional file 2: Table S3 
and Table S4). Hence, the glucose was consumed within 
15 h and ca. 80% of the xylose was consumed within 65 h 
(Additional file 1: Figure S3). However, despite the in- 
creased nutrient availability no significant differences were 
observed between the control strain and the FRD1 -over- 
expressing strains TMB3456 and TMB3458 (Additional 
file 2: Table S4). 



Over-expression of mitochondrial FR has a detrimental 
effect on xylose fermentation 

Although the change from a defined medium to a com- 
plex medium led to a general increase in the growth rate 
on xylose, the over-expression of FRD1, alone or in com- 
bination with EROl, had little effect on the cellular 
physiology. One explanation could be that the previously 
observed accumulation of fumarate [31] does not occur 
in the cytosol, but is restricted to the mitochondria. To 
investigate this possibility we over-expressed the OSM1 
gene, alone and in combination with EROl, and evalu- 
ated the fermentation performance in both defined and 
complex media with 20 g/L glucose and 50 g/L xylose. 

In defined medium, TMB3457 (OSM1) and TMB3459 
(OSM1 EROl) also consumed the glucose within 20 h 
and ca. 60% of the xylose within 90 h (Figure 4C and E). 
During glucose fermentation these strains behaved simi- 
lar to the other strains and gave nearly the same product 
yields as the control strain (Table 2). When xylose be- 
came the sole carbon source the physiology changed 
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Time (h) 

Figure 4 Fermentation of a glucose/xylose mix in defined medium. A mix of 20 g/L glucose and 50 g/L xylose was fermented in defined 
medium (2X YNB) by strains A) TMB3455 (control), B) TMB3456 (FRD1), C) TMB3457 (0SM1), D) TMB3458 {FRD1 EROl) and E) TMB3459 (OSM1 ER01). 
Symbols: diamonds, glucose; squares, xylose; plus, biomass; triangles, ethanol; circles, xylitol; stars, glycerol. Figures illustrate one representative 
experiment out of two biological replicates. 



drastically and the two strains behaved differently. The 
maximum growth rate of strain TMB3457 (OSM1) was 
47% lower (P = 0.03) compared with the control strain 
(Table 2) and the biomass-specific yields of acetate and 
ethanol increased by 130% (P = 0.04) and 69% (P = 0.05), 
respectively (Table 2). As a result a smaller fraction of 
the carbon was distributed to biomass and glycerol in 
favour of acetate and ethanol (Table 2). Strain TMB3459 
(OSM1 EROl) did not display any difference in growth 
rate on xylose compared with the control strain (Table 2). 



The acetate yield per biomass also increased by 44% (P = 
0.008) whereas the glycerol and xylitol yields decreased by 
39% CP = 0.08) and 20% (P = 0.002) compared with the 
control strain (Table 2). Hence, this strain diverted the 
carbon away from glycerol, xylitol and ethanol for a higher 
yield of acetate and maintained biomass yield (Table 2). 

In complex medium strains TMB3457 (OSM1) and 
TMB3459 (OSM1 EROl) consumed glucose within 15 h 
but only 50-60% of xylose was consumed within 65 h 
(Additional file 1: Figure S3). The performance of the 
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Table 2 Maximum specific growth rates and product yields in anaerobic batch fermentation of 20 g/L glucose and 



50 g/L xylose 


in defined medium (2X YNB) using strains TMB3455, TMB3457 and TMB3459 








VJIUV.UOC yji icijc 






Ay iujc yji icijc 






TMB3455 
(control) 


TMB3457 
(OSM1) 


TMB3459 
(OSM1 EROl) 


TMB3455 
(control) 


TMB3457 
(OSM1) 


(OSM1 EROl) 


Umax Ufa) 


0.29 ±0.01 2 


0.30 ±0.01 9 


0.32 ± 0.044 


0.029 ± 0.002 


0.01 8 ±0.001 


0.028 ± 0.003 


Yields per substrate (mol/mol sugar) 












Yxylt/s 




- 


- 


0.24 ± 0.03 


0.25 ±0.01 


0.20 ± 0.02 


Yglyc/s 


0.22 ± 0.02 


0.20 ± 0.02 


0.20 ± 0.05 


0.020 ± 0.000 


0.01 6 ±0.002 


0.01 2 ±0.001 


Yac/s 


0.020 ±0.001 


0.01 2 ±0.000 


0.01 7 ±0.005 


0.01 9 ±0.001 


0.031 ±0.003 


0.029 ± 0.003 


Yetoh/s 


1.87 ±0.14 


1.96 ±0.21 


1 .96 ± 0.37 


1 .03 ± 0.00 


1.24 ±0.1 5 


0.87 ± 0.00 


Ysucc/s 


0.003 ± 0.000 


0.001 ±0.000 


0.002 ±0.001 


0.009 ± 0.000 


0.007 ± 0.002 


0.006 ± 0.000 


Y x/S (g/mol) 


10.4 ±0.5 


10.5 ±0.2 


1 1 .9 ± 0.9 


3.8 ± 0.4 


2.7 ± 0.6 


3.9 ± 0.4 


Yields per biomass (mmol/g CDW) 












Yxylt/x 








63.1 ±0.3 


92.9 ±25.8 


50.3 ± 0.8 


Yglyc/x 


20.2 ± 0.8 


18.3 ±0.5 


16.6 ±0.2 


5.2 ± 0.6 


5.8 ±0.7 


3.2 ± 0.6 


V 

1 ac/x 


2.1 ±0.1 


1.8 ±0.1 


1 .8 ± 0.0 


5.1 ±0.3 


1 1 .7 ± 1 .8 


7.3 ± 0.0 


Y etoh/x 


153 + 8 


153±11 


147±7 


271 ±30 


458 ±53 


223 ± 23 


Ysucc/x 


0.1 9 ±0.04 


0.28 ±0.1 8 


0.22 ± 0.04 


2.5 ± 0.4 


2.4 ± 0.3 


1 .6 ± 0.22 



Values are given as mean ± standard deviation of two independent experiments. Data for strains TMB3456 and TMB3458 is reported in Additional file 2: Table S3. 



strains did not differ significantly from the control strain 
during the glucose phase (Table 3). With xylose as the 
carbon source, on the other hand, both strains grew 
markedly slower than the control strain. The growth 
rates of TMB3457 (OSM1) and TMB3459 (OSM1 EROl) 
on xylose were reduced by 25% (P = 0.09) and 51% 



(P = 0.12), respectively (Table 3). Both strains also in- 
creased the biomass-specific yield of acetate by 73% 
(P = 0.003) and 157% (P = 0.03), respectively (Table 3). 
Although both strains diverted more carbon toward acet- 
ate, different by-product yields were reduced by each 
strain. TMB3457 (OSM1) gave a reduced substrate-specific 



Table 3 Maximum specific growth rates and product yields in anaerobic batch fermentation of 20 g/L glucose and 



50 g/L xylose in complex medium (10 g/L yeast extract) using strains TMB3455, TMB3457 and TMB3459 






Glucose phase 






Xylose phase 






TMB3455 


TMB3457 


TMB3459 


TMB3455 


TMB3457 


TMB3459 




(control) 


(OSM1) 


(OSM1 EROl) 


(control) 


(OSM1) 


(OSM1 EROl) 




0.37 ± 0.003 


0.38 ±0.001 


0.38 ± 0.003 


0.033 ± 0.004 


0.024 ± 0.000 


0.01 6 ±0.008 


Yields per substrate (mol/mol sugar) 












Yxylt/s 








0.14 ±0.01 


0.1 6 ±0.01 


0.1 3 ±0.01 


Yglyc/s 


0.13 ±0.01 


0.11 ±0.01 


0.11 ±0.01 


0.01 3 ±0.001 


0.009 ±0.001 


0.009 ± 0.002 


YqcA 


0.030 ± 0.006 


0.022 ± 0.000 


0.025 ± 0.003 


0.022 ± 0.003 


0.037 ± 0.005 


0.037 ± 0.001 


Yetoh/s 


2.7 ± 0.6 


2.6 ± 0.2 


2.6 ± 0.2 


1.4 ±0.1 


1 .4 ± 0.0 


1.1 ±0.0 


Ysucc/s 


0.005 ± 0.004 


0.003 ±0.001 


0.003 ±0.001 


0.003 ± 0.001 


0.006 ± 0.000 


0.005 ± 0.001 


Y x/S (g/mol) 


17.5 ± 1.4 


18.4 ±0.9 


19.5 ±0.6 


6.4 ± 0.8 


6.3 ± 0.7 


4.2 ± 0.8 


Yields per 


biomass (mmol/g CDW) 












Y X ylt/x 








22.3 ± 0.8 


25.5 ±1.2 


31.0 ±8.7 


Yglyc/x 


7.3 ± 0.5 


6.4 ±0.2 


6.4 ± 0.0 


2.1 ±0.0 


1 .5 ± 0.4 


2.3 ± 0.9 


Yqc/x 


1 .3 ± 0.4 


1.3 ±0.1 


1.5 ±0.1 


3.4 ± 0.0 


5.9 ±0.2 


8.8 ± 1 .4 


Y etoh/x 


94.2 ± 5.0 


91.3 ±5.9 


93.0 ±3.5 


222 ± 1 1 


228 ±21 


272 ±62 


V 

1 succexxx 


0.12 ±0.07 


0.1 7 ±0.08 


0.21 ±0.02 


0.54 ± 0.03 


0.91 ±0.03 


1 .2 ± 0.5 



Values are given as mean ± standard deviation of two independent experiments. Data for strains TMB3456 and TMB3458 is reported in Additional file 2: Table S4. 
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yield of glycerol by 31% (P = 0W) whereas TMB3459 
(OSM1 EROl) gave reduced yields of ethanol and bio- 
mass by 21% (P = 0.06) and 34% CP = 0.11), respectively 
(Table 3). 

These results indicate that the reactions taking place 
in the mitochondria have a significant role in xylose me- 
tabolism and that they are more sensitive to changes in 
metabolite levels (leading to altered regulatory signals, 
degree of allosteric inhibition of enzymes etc.) on xylose 
compared to glucose. Hence, the over-expression of 
OSM1 caused physiological responses that were not ob- 
served in strains over-expressing FRD1. In strain 
TMB3457 (OSM1) these responses were: i) a 20-40% re- 
duction in glycerol yield, ii) a 63-68% increase in acetate 
yield and iii) a 27-38% reduction in growth rate during 
xylose fermentation (Table 2 and Table 3). A possible ex- 
planation for this phenotype could be that Osmlp inter- 
acts with Gut2p and Achlp which are both located 
specifically inside the mitochondria (Additional file 1: 
Figure S4). GUT2 encodes a FAD-dependent glycerol 3- 
phosphate dehydrogenase located in the mitochondrial 
inner membrane [41] and expression of the gene is re- 
pressed in the presence of glucose and induced during 
growth on non-fermentable carbon sources [42]. An in- 
duction of the GUT2 gene has also been shown to occur 
during xylose fermentation with transcript levels more 
than 7-fold higher compared to glucose [27] (Additional 
file 2: Table S5). Recently, the protein encoded by ACH1 
was shown to function as a Coenzyme A (CoA) transfer- 
ase rather than an acetyl-CoA hydrolase [43]. Hence, this 
enzyme catalyses the following reaction: acetyl- CoA + suc- 
cinate — > succinyl-CoA + acetate (A r G'° = -8.1 kj/mol) and 
is also induced nearly 5 -fold during xylose fermentation 
compared to glucose [27] (Additional file 2: Table S5). If 
both Gut2p and Achlp are active on xylose, the pathway 
formed together with Osmlp could explain the diversion 
of carbon from glycerol to acetate as well as the reduced 
growth rate due to a reduced acetyl-CoA availability 
(Additional file 1: Figure S4). However, this hypothesis 
needs to be validated through further investigations. 

Frdlp and Osmlp do not influence regeneration of 
cytosolic NAD 

An early study on the role of FR showed that the 
addition of oxidised methylene blue or phenazine meth- 
osulfate, which chemically oxidises NADH to NAD, to 
cultures of a frdlA osmlA double mutant rescued the 
non-growing phenotype under anaerobic conditions 
[17]. It was thus proposed that the two FR enzymes are 
involved in the oxidation of excess NADH formed in the 
cytosol and the mitochondria under anaerobic condi- 
tions. The formation of xylitol during xylose fermenta- 
tion has been ascribed to an inability of the cells to 
regenerate the NAD needed for the conversion of xylitol 



to xylulose by XDH [44] . The addition of acetoin, a sub- 
strate which can be reduced by yeast to 2,3-butanediol 
using NADH, to cultures fermenting xylose reduces the 
xylitol yield significantly [45,46]. If the FR enzymes par- 
ticipate in the conversion of excess NADH a reduction 
of the xylitol yield during xylose fermentation would be 
expected when the FR activity is increased. However, the 
over-expression of FRD1 or OSM1 did not affect the 
amount of xylitol produced from xylose (Table 2, Table 3, 
Additional file 2: Table S3 and Table S4). This is in 
agreement with other studies showing that addition of 
acetoin to a culture of a frdlA osmlA double mutant 
does not rescue the non-growing phenotype under an- 
aerobic conditions [11]. In addition, the proposed mech- 
anism of NADH oxidation requires the presence of an 
enzyme capable of transferring electrons from NADH to 
FAD to generate the reduced cofactor used as substrate 
by fumarate reductase [11,17]. As of yet, no such en- 
zyme has been identified in S. cerevisiae that is ex- 
pressed under anaerobic conditions. There are, however, 
FR enzymes that use NADH as co-factor instead of 
FADH 2 . The NADH-dependent FR enzyme from Trypa- 
nosoma brucei has recently been expressed in a xylose- 
utilizing S. cerevisiae strain which resulted in increased 
ethanol yield and decreased xylitol yield when ferment- 
ing xylose [47]. Based on these results we conclude that 
the FR enzymes from S. cerevisiae do not influence the 
regeneration of cytosolic NAD. 

Other potential limitations in the secretory pathway in 
S. cerevisiae during xylose fermentation 

The results obtained in this study indicate that the activ- 
ities of FR and Erolp enzymes are not limiting during xy- 
lose fermentation. Measurements of the FR activity in 
strain TMB3455 when growing on glucose and when me- 
tabolizing xylose indeed showed equal activities around 
3 mU/mg protein in both conditions (Additional file 1: 
Figure S5). This is contradicting previous studies showing 
that the activity and expression of FRD1 and OSM1 de- 
crease upon entry into stationary phase [11,39] and during 
carbon starvation [48]. It is, however, possible that the FR 
activity is maintained high on xylose in response to cellu- 
lar stress. The up-stream regions of the FRD1 and OSM1 
genes contain binding sites for stress -related transcription 
factors Yaplp and, in the case of FRD1, Msn2p/Msn4p 
[49]. Due to the involvement of Yaplp, an oxidative stress 
would be most likely in this case [50]. This possibility is 
supported by much higher intracellular concentrations of 
both reduced and oxidized glutathione [31] and higher 
transcript levels of MSN4 and YAP1 [27] in recombinant 
S. cerevisiae strains during xylose fermentation compared 
to glucose fermentation (Additional file 2: Table S5). 

The up-regulation of RIB1 (Figure 3) could be an indi- 
cation that the need for flavin-containing complexes 
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and/or the free co-factor is increased on xylose. Erolp 
has been shown to use free FAD as electron acceptor 
in vitro under anaerobic conditions [51] which is in line 
with the ability of the FR enzymes to use free FADH 2 as 
cofactor and their essential nature under such conditions. 
As S. cerevisiae possesses one gene encoding a mitochon- 
drial FAD transporter (FLX1) [52] and three genes encod- 
ing putative transporters associated with the ER (FLC1-3) 
[53], FAD-exchange between compartments should in- 
deed be possible. This highlights free FAD as a potentially 
crucial metabolite for sustaining anaerobic growth. 

Pdilp is, together with Erolp, essential for disulphide 
bond formation in the ER [54]. Although transcriptional 
analysis did not show any difference in expression of 
PDI1 between glucose and xylose fermentation [27] 
(Additional file 2: Table S5), proteome analysis of a mu- 
tant strain with improved xylose-fermenting capability 
identified Pdilp as significantly more abundant in the 
improved strain compared with the wild-type [55]. This 
strongly indicates that there indeed is a need for a higher 
capacity in the protein folding mechanism for efficient xy- 
lose fermentation. Whether this capacity is for regulation 
of disulphide bond formation [56], reduction of non- 
native disulphide bonds [57] or targeting unfolded glyco- 
proteins for degradation [58] remains to be investigated. 

Conclusions 

Over-expression of FRD1, encoding the cytosolic FR, did 
not improve the growth rate on xylose, even when the 
strain was provided with all amino acids and an additional 
copy of the EROl gene controlled by a strong promoter. 
Over- expression of the OSM1 gene, encoding the mito- 
chondrial FR, even had a detrimental effect and caused a 
reduced growth rate. Hence, we conclude that increasing 
the activities of the FR enzymes and Erolp is not sufficient 
to increase the anaerobic growth on xylose. However, our 
transcription analysis of UPR related genes showed a wide 
down-regulation of the secretory pathway during xylose 
fermentation indicating that additional components of the 
protein folding mechanism may be limited. Two immedi- 
ate factors should notably be investigated further: i) the 
actual form of HAC1 inside the cells during xylose fer- 
mentation to determine the regulatory component and ii) 
the activity of Pdilp which is acting together with Erolp 
in the formation of disulphide bridges. Furthermore, un- 
derstanding the requirement of free FAD for anaerobic 
growth could prove to be essential for the development of 
efficient xylose-utilizing S. cerevisiae strains. 

Additional files 



Additional file 1: Supplementary Figures S1-S5. Description: FR 
activity measurements, sensitivity results, fermentation profiles. 



Additional file 2: Supplementary Tables S1-S5. Description: Genes 
included in the cluster analysis, maximum specific growth rates and 
product yields using strains TMB3455, TMB3456 and TMB3458. 
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